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ABSTRACT 

We have modeled a small sample of Seyfert galaxies that were previously identified 
as having simple X-ray spectra with little intrinsic absorption. The sources in this 
sample all contain moderately broad components of Fe K-shell emission and are ideal 
candidates for testing the applicability of a Compton-thick accretion-disk wind model 
to AGN emission components. Viewing angles through the wind allow the observer 
to see the absorption signature of the gas, whereas face-on viewing angles allow the 
observer to see the scattered light from the wind. We find that the Fe K emission line 
profiles are well described with a model of a Compton-thick accretion-disk wind of solar 
abundances, arising tens to hundred of gravitational radii from the central black hole. 
Further, the fits require a neutral component of Fe Ka emission that is too narrow to 
arise from the inner part of the wind, and likely comes from a more distant reprocessing 
region. Our study demonstrates that a Compton-thick wind can have a profound effect 
on the observed X-ray spectrum of an AGN, even when the system is not viewed through 
the flow. 

Subject headings: galaxies: active - X-rays: galaxies - accretion, accretion disks 

1. Introduction 

Active galactic nuclei (AGN) have rich spectra with a multitude of absorption and emission 
signatures, which offer insight into the structure of the nuclear environs. Current X-ray observa- 
tories, Suzaku, Chandr a, and XMM-Newton, can resolve H- and He-like species of C, N, O, Ne, 



Mg, Al, Si, and S (e.g.. iBrandt Sz Kaspil l2002), Absorption lines from these species generally show 
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1000 km s^. iBlustin et al.ll2005h . indicative of an origin m 



blueshifts (typically in the range 100 
outflowing gas. 

In addition, features from a range of ionization states of Fe have been resolved including the 
neutral Fe Ka emission line at ~6.4 keV; the H-like and He-like species at higher energies, which 
are often seen in e mission and absorp tion; a broad unresolved transition array from the M-shell just 
below 1 keV (e.g. iBehar et al.ll200ll ); and Fe L transitio ns, which are usually seen in absorption 
and sometimes weakly in emission (e.g. iKaspi et al.ll2002l ). Observations of a strong Fe L emission 
signature ar e sometimes evident in AGN spectra but generally attrib uted to starburst activity, e.g. 
NGC 1365 (|Guainazzi et aljboosl ) and Mrk 573 (jBianchi et aliboioh . 



In principle, the profile of Fe Ka emission can be used to constrain conditions in the repro- 
cessing gas. The line is sufficiently strong that useful constraints may be obtained from currently- 
available data. However, interpretation of the line profile has been challenging because of the 
difficulty in separating line and continuum photons in X-ray data. So me early work sugge sted 
the line to be extr emely broad, owing to strong relativistic effects (e.g. iRevnolds et al.lll995l and 
Fabian et al.l 120021 ) . Such an interpretation required the line be produced within a few r g of the 
black hole. An A SCA sample study suggested relativistically-blurred lines exi st in ~ 75% of loca l 



Seyfert galaxies (INandra et al.1 Il997l ) , although later work using XMM data (INandra et al.1 120071 ) 



concluded that only ~ 40% of local Seyferts showed evidence for lines orig inating within 50r 9 . Sub 



sequently, analysis of a flux- limited sample (jde La Calle Perez et al.ll2010l ). reduced the estimate of 
extremely broad Fe Ka emission to ~ 30% of sources. However, strong evidence for high column 
X-ray absorbers in the local AGN population has motivated several teams to reconsider the true 
line profile and origins of the Fe K-shell emission, and it may be that all claimed relativistically- 
blur red lines are actua lly a manifestation of absorption and scattering through circumnuclear gas 
(e.g. iMiller et alJbood ). 



In a study of NGC 3783, iReeves et al.l (J200J) found that modeling all of the absorbing zones 
accounted for much of the curvature originally attributed to the broadness of the Fe K emission 
line. In NGC 4051, detailed modeling of the absorber complex l ed to a limit FWH M < 5540 km 
s _1 for the moderately broad component of Fe K emission (e.g. iTurner et al.ll2010i ). In contrast, 
when partial- covering absorption was not accounted for in model fits, the mean FWHM line width 
suggested by INandra et al.l (|2007l ) was ~ 37000 km s _1 for a sample of Seyfert 1 galaxies. When 
the absorbing complex is fully accounted for there remains no requirement for a relativistically - 
blurred component of Fe Ka emission in any source, not even MCG-6-30-15 (jMiller et al.l 12008) . 
However, a very modest broad component is sometimes evident that does not require the invocation 
of relativistic effects, and this is the subject of our paper. 



To explain th e grating spectra of the best studied sources (e.g. MCG-6-30-15, ILee et al 



2001 



Miller et al.l 120081 ) . multiple zones of ionized, absorbing gas are required. Of particular interest is 



cm 



Pounds et al. 
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Risaliti et al. 
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20051 : iMiller et al.l 120071 ; iTurner et al.l 120081 ) . Such high-column absorbing system s appear to be 



comm on in radio-quiet AGN, being detected in more than 40% of sources studied bv lTombesi et al 
(201CJ) , and typically showing significant outflow velocities from a few thousand km s _1 upwards 



to ~0.3c. 



The strong evidence for Compton-thick outflows in type 1 AGN motivat ed development of a 



mult i -dimen sional radiative-transfer realization of a Compton-thick wind (e.g. ISim 



2008, l2010al ). Such winds are expected in AGN having a high Eddington rat io (IKing fc Pounds 



20031 ) and these would likely be launched fairly close to the black hole (e.g. iProga et al 



2005; 



Sim et al 



200' 



0) 



Consideration of high wind densities led to the prediction of strong wind signatures, including deep 
absorption lines and a moderately broadened Fe K-shell emission profile. 



Sim et al.l (J2008J, l2010al ) have calculated wind spectra assuming a launched radius beyond 30 
gravitational radii, such that relativistic blurring does not significantly affect the line width. These 
simulations show a moderately broadened line profile having FWHM ~ 1 keV owing to the combined 
effects of Compton scattering and the large rotational/outflowing velocities in the wind. Electron 
scattering in the flow combined with photoelectr ic absorption contributes to a slightly red-skewed 
wing to the observed Fe K emission profile (c.f. lLaurent &; Titarchuki 120071 ) . The Fe Ka line is a 
strong signature in the X-r ay spectrum of the wind, with equivalent widths in excess of ~ 200 eV 
for some wind parameters (|Sim et al.ll2008l ). 



Thus, a Compton-thick reprocessing wind simultaneously produces both blue-shifted absorp- 
tion lines and red-skewed emission lines in the Fe K regime for AGN, the relative strength of these 
in the observed spectrum depends upon the orientation of the flow with respect to the observer's 
line-of-sight. 

In this paper, we focus on the interpretation of observable features in the Fe K band in the 
context of a Compton-thick wind. Although other spectral lines are expected to form in the flow, 
identification and isolation of features in the Fe K energy band is more secure than in the soft 
X-ray regime, since only K shell transitions of heavy ions are expected to be present around 6 - 
7 keV. Although Fe K-shell absorption features are excellent diagnostics for a highly ionized flow, 
the associated Fe K emission arising from line scattering or reco mbination i n the outflow is also 
important and has only been assessed for two AGN to date (see ISim et al.l l2010al ) . Specifically, 
case studies of Mrk 766 and PG1211+143 have shown that absorption and emissi on features in the 



2010a 



Fe-K regime can successfully be explained using a single Compton-thick flow (see ISim et al. 
and references therein). 

The best constraints on the emission line profile can be obtained when there are no confusing 
absorption features present; these are diffi cult to resolve s pectr oscopically using current data. With 
this in mind, the small sample identified bv lPatrick et al.l (|201ll ) is of particular interest; that sample 
comprises objects showing little evidence of intrinsic X-ray absorption. Unobscured objects present 
a relatively clear view of the Fe K region enabling a fit to the broadband spectrum and a basic 
understanding of the Fe K emission features that is unencumbered by consideration of separating 



- 4 - 



the signature of a complex X-ray absorber. 



Patrick et al.l (|201ll ) found a moderately broad component of Fe Ka to exist in all of their 



sample sources. Here we test whether these moderately broad Fe-K emission, together with limits 
on absorption signatures, components are consistent with an origin arising from reflection off a 
Compton-thick wind. We begin in Section 2 by outlining the parameters of the wind model. In 
Section 3, we outline the observations and data reduction; and in Section 4, we detail the spectral 
fitting process and results. Section 5 discusses the implications of our results. We draw conclusions 
in Section 6. 



2. The Compton-thick Wind Model 



To model the Compto n-thick wind, we used a model grid based on the Compton-thick wind 
model of ISim et al.l (|2010al ). We present a general overview of the Compton-thick wind model and 
refer the reader to the aforementioned paper for specific details. 



Sim et al.l (|2010al ) presented the results of Monte Carlo radiative transfer calculations, sim- 



ulating a Compton-thick disk wind. The wind was parameterized as a steady-state, continuous, 
bi-conical outflow. The density in the flow varied smoothly as a function of position and was de- 
termined by the geometry of the flow, the total mass-loss rat e (M,„) and the relative mass-loading 
of streamlines in the flow. Hydrodynamical simulations (e.g. iProga &: KaHmanll2004l ) showed that 
disk wind s were e xpected to have very complex, non-spherically symmetric geometries. However, 
Sim et al.l (|2010al ) invoked a simple velocity parameterization because it allowed for a relatively 
simple geometry. The model does not claim to offer a completely realistic description of the flow, 
but approximates the key physical properties of the wind and makes it possible to unde rstand the 



influence of these parameters on the observed spectrum, by computing model grids (see ISim et al 



2008 



2010al ). Compari son between this version of the model and points from a more complex, 



hydrodynamical model (|Sim et al.ll2010bl ) suggest that the velocity prescription can yield predic- 
tions for the wind that are close to those obtained with a much more time-consuming, intensive 
calculation. 



The geometry adopted a modification of the "displaced dipole" model described in lKnigge et al 



(j 19951 ). Three parameters defined the extent and opening angle of the wind: was the distance below 
the origin to which all the stream lines in the wind converge; R m in was distance from the origin 
to the innermost edge of the wind in the xy-plane in units of gravitational radius, r 9 ; and R max 
was the distanc e from the origin to the outermost edge of the wind in the xy-plane (see Figure 1 in 
Sim et al.ll2008l ). The disk was assumed to lie in the xy-plane and the wind was symmetric under 
rotation about the z-axis and under reflection in the xy-plane. 



Sim et al.l (|2010al ) invoked a p rescription th a t spec ifies the velocity at every point in the wind 



following the parameterization of iKnigge et al.l (|1995l ) . It was assumed that parcels of matter 
conserve specific angular momentum about the z-axis as they flowed outwards with the angular 
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momentum at the base of a streamline set at the Keplerian value for the radius at which this 
rose from the xy-plane. Thus, the rotational velocity was determined by the specification of wind 
geometry and the mass of the central black hole. The poloidial velocity equation was 

vi=v + (voo - v )(l - - — ■ y 

Jx v + I 

where = fv\J^$^ was a user specified multiple, f v , of the escape velocity from the base of 
the streamline in the xy-plane, Mbh was the black hole mass, Vq was the initial velocity, I was the 
distance along the outflow streamline, R v was an acceleration length parameter, and (3 set the rate 
of acceleration. Varying the /?, R„ and f„ parameters allowed us to explore the influence of different 
rates of wind acceleration in the model. Ultimately, this related to the physical process responsible 
for driving the flow. 



As discussed bv lSim et al.l (120081'). th e Monte Carlo simulations relied on discrete, indivisible 



packets of radiation based on iLucyl ([2002). In order to calculate the radiative transfer, a 100 x 100 
wind grid was produced where the wind properties was discretized in terms of the parameters R v 
and I. This grid was transformed into a 3D 100 x 100 x 100 Cartesian grid, where those grid points 
inside the wind were assigned density, ionization/excitation states and radiation field properties 
of the point closest to the wind grid point. Those grid points outside the wind remained empty. 
These discrete packets ranged from 0.1 keV to 511 keV with an assumed power-law distribution 
with the photon index (T) as an input parameter. The X-ray source was assumed to be a spherical 
shell with radius 6 r g , so that the X-ray emitting source had a non-zero angular size. The packets 
propagated until they reached the outer boundary of the computational domain. They may have 
interacted with the wind via Compton scattering, bound-bound, bound-free (including photoelectric 
absorption by inner shell elec trons) and free-free processes. Following interactions, the Macro Atom 



method of iLucvl ([2002I . 120031 ) was used to simulate subsequent fluorescent or thermal re-emission 
by the outflow. During each Monte Carlo simulation, the wind properties (e.g. ionization state, 
temperature) were held fixed. Once all packets propagated through the outflow, the packet histories 
were used to compute photo-ionization and heating rates in each grid cell. These were then used 
to solve the equation of ionization and thermal equilibrium and determine improved estimates for 
the ionization state and temperature. These estimates were then used as input to a new Monte 
Carlo simulation. Thus, via a sequence of Monte Carlo simulations, we iteratively determined the 
conditions in the outflow. The total output spectra contained two components: a direct photon 
spectrum, containing photons that had no interaction with the wind, and a scattered/reprocessed 
photon spectrum, containing photons that interacted with the wind. The output spectrum was not 
calculated until the last iteration of the simulation. 



2.1. Construction of synthetic wind spectra for spectral fitting 



To generate a set of synthetic wind spectra to use for fitting the data, we performed radiative 
transfer calculations for a grid of models. In the grid, we varied the following parameters of the 
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model: M w , d, R m m ; Rmsi, Fe abundance and T. The values used on the grid are listed in Table[TJ 



The spectral signature from the Compton-thick wind model of lSim et al.l (|2010al ) was calculated 
over a range of parameter-space. We assumed f v = 1, which is consistent with a radiatively driven 
flow. The parameters R v and (3 depend on the physical acceleration process and the wind location; 
we assumed R v = R m ax in ah cases and (3 = 1, and these represent acceleration which occurs on 
the scale of the system. We assumed vo <C , and therefore adopted v o = . A black hole mass 
M= 10 7 Mq was utilized, with an intrinsic 2-10 keV luminosity Lx = 0.01 l>Edd, where l^Edd was 
the Eddington luminosity, and a mass-loss exponent k= -1, which was chosen because this gives a 
gentle density variation with radius making the ionization structure of the wind easier to resolve. 
The outer radius of the simulation grid, i.e. the extent of the region over which the influence of the 
wind on the spectrum was calculated, was 5 x 10 16 cm. For models with higher Fe abundance, we 
also increased the Ni abundance by the same factor. 

Given the limited spectral resolution and signal-to- noise ratio of the data, we found that it 
was not possible to determine all wind parameters simultaneously. As a result of these parameter 
degeneracies, we fixed both the inner radius and the focal point d at 32 r p ; the outer radius to ~ 



300 T g and the Fe abundance to solar abundance, via lAsplund et al.l (|2005l ). This geometry fixed 
the opening angle to ~ 49°. We further discuss the parameter degeneracies in Section 5. Although 
a black hole mass is specified, the simulated wind is applicable to other black hole mass values as all 
parameters scale with black hole mass. This particular black hole mass was chosen as a convenient 
scaling factor. 

For each model on the grid, we computed synthetic spectra for 20 observer orientations. We 
therefore generated a total of 4320 synthetic spectra. 



3. Observations 

Suzaku has four X-ray telescopes, each containing a silicon CCD within its focal plane forming 
the X-ray Imaging Spectrometers (XIS) suite. XIS0, XIS2 and XIS3 are front-illuminated (FI), 
providing data over a usable range of 0.6-10.0 keV with an energy resolution FWHM ~ 130 keV at 
6.0 keV. In November 2006, a charge leak was discovered in XIS2, making XIS0 and XIS3 the only 
operational FI chips. XIS1 is back-illuminated. The back-illuminated configuration extends the 
soft band to ~0.2 keV; however, it also results in a lower effective area and higher background rate 
in the Fe K regime, compared to the FI chips. Consequently, XIS1 is excluded from our spectral 
analysis. Suzaku also carries the Hard X-ray Detector (HXD) that contains a silicon PIN diode 
detector covering a range of 10.0-100 keV with a usable energy range of 15.0-70.0 keV. 

The Suzaku observations were made 2006 August 28 (OBSID 701031010), 2007 April 1 (OBSID 
702014010), 2007 June 7 and 2010 May 19 (OBSID 702043010, 705063010), 2007 December 9 
(OBSID 702122010), 2008 June 24 (OBSID 703028010), and 2008 August 28 (OBSID 703060010) for 
Mrk 335, Ark 120, Fairall 9, SWIFT J2127.4+5654, NGC 7469, and MCG-02- 14-009, respectively. 
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An observation log, observation exposure times and source count rates are detailed in Table 2 of 



Patrick et al.1 (|201lh 



The data were reduced using HEAsoft v. 6. 10. The XIS cleaned event files were screened in 
XSELECT to exclude data during passage through the South Atlantic Anomaly and also excluding 
data starting 500 s before entry and up to 500 s after exit. In addition, we excluded data having 
an Earth elevation angle < 10° and a cut-off rigidity > 6 GeV. CCDs were in 3 x 3 and 5x5 edit 
modes, with normal clocking mode. Good events were selected, having grades 0, 2, 3, 4, and 6, 
while hot and flickering pixels were removed using the SIS CLEAN script. The spaced-row charge 
injection was utilized. XIS spectra were extracted from circular regions of 3.0' radius centered on 
the source, while the background was extracted from a region of the same size offset from the source 
and from the corners of the chip that register calibration data. 

Data from the front-illuminated CCDs XIS0 and 3 were used for all observations, and data 
from the two chips were were co-added using addascaspec. For Mrk 335, data were also available 
from XIS2, and these were combined with XIS0,3 in that case. For SWIFT J2127. 4+5654, only 
data from the XIS3 chip were used, since the XIS0 chip was not functioning properly at the time 
of observation. 

The cleaned PIN data was reduced utilizing the ftool hxdpinxbpi. This tool calculates good 
time intervals (GTIs) of the non X-ray instrumental background (NXB, using model 'D' released 
2008 June 170 ) data that overlap with the source data and extracts both the source spectrum 
and NXB background spectrum through that common GTI. A simulated cosmic X-ray background 
(CXB) spectrum is calculated and combined with the NXB spectrum to produce a total PIN 
background spectrum for the observation. A dead-time correction ~ 4-5% is applied to the source 
spectrum using the ftool hxddtcor. Finally, the PIN data have a known la systematic uncertainty 
of 1.39c@, which was applied to the PIN data in grppha. 



4. Spectral Fitting Results 

Spectral analysis was performed using XSPEC v 12.5. To fit these data, the 4320 syn- 
thetic spectra previously mentioned were used to create a multiplicative model table (mtable). To 
adjust for instrumental calibrations during spectral fitting, the PIN model was scaled by an energy- 
independent constant as appropriate for the HXD and XIS nominal aim points, respectively. The 
constant value was 1.18 for all observations except Fairall 9 (OBSID 702043010), for which a con- 
stant scalar of 1.16 was appropriate. Prior to fitting, the XIS data were binned to the half- width 
half maximum resolution of the instrument, where as the PIN data were binned to Alogio Energy= 
0.1 resolution. In the fitting process, the \ 2 statistic was minimized and then errors were calculated 



1 http://www.astro. isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2007-01.pdf 
2 htt p: / /heasarc. nasa. gov / docs /suzaku / analysis / watchout . html 
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for the 90% confidence interval for each parameter. 

First, we fit the 2-10 keV band, to investigate whether the wind model can account for the Fe 
K-shell spectral signature in the sample sources. We note that the two observations of Fairall 9 were 
fit simultaneously to constrain parameters that should not differ between epoch, such as viewing 
angle and Galactic column. The spectral model was constructed using a powerlaw continuum of 
free normalization and photon index T. The continuum form was modified by the absorption and 
remission from the Compton-thick wind and the illuminating spectrum for the wind was assumed 
to be the power-law continuum. As stated previously, the degeneracies in the model have led us to 
allow for only three free parameters in the wind model: T, M w and 9. 



Patrick et al.l (120111 ) previously identified several components of Fe Ka emission in the sample 
sources: most objects show clear moderately broad and narrow components of Fe Ka emission. To 
account for the core of the Fe emission line, which is generally too narrow to be produced in the 
wind, we included a narrow (<r= 10 eV) Gaussian model component. All models also includ ed a 



Galactic component of absorption parameterized using the tbabs model (jWilms et al.l 12000) and 



using values from the we ighted average Njj in the Dickey & Lockman HI in the Galaxy survey 



(jDickev Lockmanlll990r ) 



Our initial model consisting of the wind table, a narrow Gaussian component and tbabs 
provided a good fit to the overall spectral curvature and, in particular, to the moderately broad Fe 
K emission component of the sample sources. However, the model systemically underestimated the 
depth of the Fe K absorption edge. Adding the PIN data to the fit and extrapolating the model 
up to 50 keV revealed that this construction also underestimated the PIN flux. Therefore, we 
considered whether adding a cold, neutral reflector would improve the fit. The addition of a cold, 
neutral reflector (f ixed at R= -1) component to the total model (parameterized using pexmon, 



Nandra et al.ll2007l ) produced an improvement to the fit to both the Fe K regime and the hard-band 
data (Figures Q] and [2]). pexmon includes the Fe Ka, Fe K/3 and Ni Ka line emission expected 
from a reflecting slab, and no additional Gaussian component was required to explain the narrow 
core of Fe K emission. As a result, we removed the Gaussian component from the total model. Our 
final model construction for fits above 2 keV consisted of the wind model, TBABS, and pexmon. 

With the addition of the cold, neutral reflector, the fits improved with A% 2 > 34 in 5/6 
sources compared to fitting the 2-50 keV data without this component. The reflector characterized 
by pexmon constituted ~ 10 - 30 % of the PIN flux, whereas the transmitted and scattered 
components of the wind constituted ~ 50 - 85 % of the PIN flux. The component contributions 
were measured against the total model predicted 15-50 keV band count rate. When statistically 
evaluating the merit of adding another component with the F-test, we found that the addition of the 
neutral reflector yielded an F-test probability of < 10~ 5 in 5/6 objects. This probability determines 
whether the null hypothesis can be rejected. Such a low F-test probability offers strong evidence 
for the addition of the extra component. We note that Ark 120 and Fairall 9 OBSID 705063010 
required an additional component of Compton-thick gas (2 x 10 24 cm -2 ) in the line-of-sight to fit 
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the PIN data, with a 12% and 31% contribution to the PIN-band flux and a covering fraction of 
8% and 18%, respectively. 

Despite the good statistical errors found for the wind model values, we note that the determina- 
tion of errors is unlikely to be better than the coarseness of the model grid. The linear interpolation 
between grid points may yield solutions that do not accurately represent the true spectral shape of 
the model in those regimes, which in turn would yield incorrect wind model parameter values. 

To test the accuracy of solutions found using linear interpolation, we compared the model 
shape for one of the fits with a fine model grid calculated near the fit solution. The wind parameter 
values using this finer grid were within the la errors of the interpolated values with a comparable 
Xv- This suggests that our coarse model grid yields reasonable wind parameter values with a linear 
interpolation. 

Because of the degeneracies in the model, the wind model parameter values must also be taken 
in context with the assumed geometry, i.e. the goal was to demonstrate that the Compton-thick 
wind model can explain the profiles of the broad components of Fe Ka emission for a reasonable 
area of parameter-space. 



Patrick et al.l (|201ll ) identified residuals at ~ 6.7 and 6.9 keV compared to their baseline model 
(Model A in their paper) . These features were attributed to narrow (cr=10 eV) ionized emission 
components from Fe xxv and Fe xxvi, superposed on the broader spectral signatures. In our 
model, the scatterin g line from the wind contributes more flux in the 6.7-7.0 keV band than the 



Patrick et al.l (|201ll ) model , giving a different baseline for any additional Gaussian line. From our 
analysis, inclusion of additional narrow Gaussian components in the model did not improve the 
fits in 4/6 objects, typically yielding A% 2 ~ 1. However, the addition of a Gaussian component to 
the total model, statistically improved the fits (A% 2 > 10) in Mrk 335 and Fairall 9. A narrow (a 
=10 eV) Fe xxv emission line at 6.67 ±0.02 keV and a narrow (<r= 10 eV) Fe xxvi emission line 
at 6.95 ± 0.05 keV were found in Mrk 335 and Fairall 9, respectively. The equivalent widths and 
line fluxes were 14±^ 4 eV and 0.23l^| x 10~ 5 photons cm" 2 s" 1 for Mrk 335 and 24±g eV and 
0.56±o;i9 x 10~ 5 photons cm" 2 s" 1 and 13±g eV and x 10~ 5 photons cm" 2 s" 1 for Fairall 

9 OBSID 702043010 and 705063010, respectively. The Mrk 335 and Fairall 9 OBSI D 702043010 



values were found to be similar to those tabulated in Model A of lPatrick et al.l (|201ll ) 



Expanding the fits to 0.6-50 keV, we found a soft excess in 4 objects: Ark 120, Fairall 9, Mrk 
335 and NGC 7469. In these objects, we parameterized the soft excess with the COMP TT model 



which represents the Comptonization of soft photons in a hot plasma above the disk (jTitarchuk 



19941 ). A soft photon temperature of 0.02 keV was used to account for the soft excess. The plasma 



tem perature and optical depth were found to be consistent with those values tabulated in Model 



A of lPatrick et al.l (|201ll ). 



Fit parameters for the wind that were obtained over 0.6-50 keV are consistent with a majority 
of those derived from the 2-10 keV band (see Table |2|). This comparison demonstrates that the 
determination of the wind parameters is robust to the approach taken with fitting, being based 
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mainly upon the profile of the moderately broad Fe K-shell emission component. From Tabled the 
average photon index based on the 0.6-50 keV fits was T = 1.96^qq} and viewing angle 9 = 59^ 
degrees. Excluding the poorly constrained mass-loss rate of Swift J2127. 4+5654, the average mass- 
loss rate was 0.89^'^ M-Edd> where M-em is the Eddington mass-loss rate. The Eddington ratios for 
this sample (Table [3]) range from 0.01-0.03, which is consistent with the assumed Lx/ ^Edd in the 
model. In the soft band of Mrk 335 and SWIFT J2127. 4+5654, we note emission and absorption 
features in the residuals of the 0.6-50 keV fits. 



Discussion 



Patrick et al.l (|201ll ) identified this sample of six type I AGN having relatively simple X-ray 
spectra, with little evidence of any intrinsic X-ray absorption. Analysis of the sample sources 
showed them to possess both moderately broad and narrow components of Fe K-shell emission. 
Modest broad components remain ing after full considera tion of the X-ray absorber complex hay e 
also been identified in NGC 5506 (iGuainazzi et alJhoiol ) and NGC 3516 (jMarkowitz et al.ll2008l ). 
Further to these, moderately broad Fe K g emission compon ents have been fo und on the variabl e 
component of the spectrum for Mrk 766 ^Miller et al.l 120071 ) and NGC 4051 (jMiller et alihoid ). 
Thus, moderately broad (FWHM < 1 keV) components of Fe K-shell emission appear quite common 
in the Seyfert population, and provide an important diagnostic of the circumnuclear gas. 

Usi ng a relativist i c disk line model to parameterize the broadness of the Fe Ka emission com- 
ponent, |Patnck_et_a]J (|201ll ) suggested that the lines originate from tens of gravitational radii, and 
that the profiles are not significantly modified by gravitational blurring. 

An alternative explanation for the moderately b road line component is that it arises from 



a Compton-thick accretion disk wind. As shown by (|Sim et al.l 120081 . I2010al ). reprocessing in a 
Compton-thick wind produces a strong broad line. In this case, rotational and outflowing velocities 
in the wind dominate the shape of the line profile, with an additional contribution from electron 
scattering. As is the case for the relativistic diskline model, the observed line shape can, in prin- 
ciple, provide information about the reprocessor, including the mass flow rate, physical state and 
geometry of the gas. When strong He- and H-like absorption signatures are detected, the wind 



(Sim 


2005; 


Sim et al. 


2008. 


2010a) 



work, comparisons between the Compton-thick wind and the observational data were limited to 
two case studies; the Fe K-shell emission and ab sorption in Mrk 766 an d PG1211+143 have been 



shown to be consistent with arising in a wind (|Sim et al.l 120081 . l2010al ). Our study extends the 



comparison of the wind model with X-ray spectral data to an additional six sources. 

Application of the wind model results in good fits to the moderately broad component of Fe 
K-shell emission in all sample sources, i.e. the breadth of the Fe K emission is consistent with the 
Compton scattering and velocity dispersion in our Compton-thick wind (Figures [TJ [2j Table [2]) of 
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solar abundance, t ested here acro s s the radial range ~ 30 — 300 r 9 . Our conclusions are qualitatively 
similar to those of iPatrick et al.l (|201ll ) in that the moderately broad Fe K emission is consistent 
with an origin tens of r 5 from the black hole. 

Since Fe K-shell emission features are formed as a result of line scattering and recombination in 
the outflow, their strength depends on t he density of ma t erial in the wind. In the model, the lines 
increase in strength with increasing M„, ([Sim et al.ll2008l . l2010al ) and so, in principle, the observed 
line strength gives a way to measure the mass loss rate in the flow. In addition, better constraints 
can be placed on M w with higher spectral resolution in the Fe K regime. 

The tabul ated range for M„, c orresponds to ~ 2 — 20 times the accretion rate estimated for th e 



sample AGN (jWoo fc Urrvi 12002 : IPetrucci et all 12004 : IVaughan et al.l 12004 : iMiniutti et al.1 12003 ). 
assuming a radiative efficiency of ~0.1. This means that our test has been limited to the case 
where outflow rates significantly exceed the accretion rate. However, such a test is reasonable, as it 
is well established fro m UV data that mass outflow rates exceed the estimated accretion rate (e.g. 
Crenshaw et al.l 120091 ) . These fits simply represent the case where most of the infalling material 
fails to reach the center of the AGN and accrete onto the black hole. 

A core component of Fe Ko emission is present at 6.4 keV in all sources. The line core is too 
narrow to be produced in the wind, but may arise via reflection from material further out. pexmon, 
which links the narrow emission line strength to the reflected continuum, better constrained the 
contributions at 6.4 keV from a cold, neutral reflector. Using this parameterization, the calculated 
relative reflection values, R/ rac , are poorly constrained, but range from 0.12 - 4.48. In addition, 
we found a preferred inclination for the distant reflector to be 85° in most objects, based on the 
joint constraints offered by the line core and hard-band flux. Such an extreme solution suggests 
that PEXMON may not be the correct parameterization for the cold, neutral reflector. However, 
Yaqoob et al.l (|2010j) have shown that a toroidal X-ray reprocessor with Compton-thick lines-of-sight 
may have Fe Ka line efficiencies that are orders of magnitude less than the maximum efficiency for 
that geometry, e.g. the Fe Ka emission line strength/ reflected continuum ratio may be different 
than those predicted by PEXMON allowing for weaker line strengths for a given hard-band flux. 
A less constrained line to reflected continuum ratio can, presumably, reduce the need for such an 
extreme inclination angle owing to having more reflection without significant production of Fe Ka. 



In addition to the narrow line core at 6.4 keV, IPatrick et al.l (|201ll ) reported the presence of 
narrow, ionized line emission from Fe xxv and Fe xxvi in some of the sample spectra. However, 
for the Compton-thick wind, the scattered Fe K line is predominantly formed as Fe xxv and xxvi, 
modified by scattering and velocity effects in the flow. Thus, our wind mode l accounted for more 
flux in the 6.5-7.0 keV regime than the baseline model of lPatrick et al.l (120111 ). and we required no 
separate model components in that energy regime, with the exception of Fe XXV emission and Fe 
xxvi emission in Mrk 335 and Fairall 9, respectively. 



There is little evidence for Fe xxv and Fe xxvi absorption in the observed spectra of the 
sample. Spectral fit results suggest intermediate viewing angles for most of these wind systems, 
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with values clustered about the mean ane le 9 = 59t 7 7 de grees. Systems at such orient ations are 
expe c ted to present a strong contribution from the scattered component of the wind (jSim et al 



2008, l2010al ). and this, in part, explains the lack of strong predicted absorption lines in the total 
spectrum. 

Further to this, as noted previously, several parameters were fixed due to degeneracies in the 
model and the values found for 9 and M w must be take n within the contex t of the assumed geometry. 
In comparison to the model parameters assumed in ISim et al.l (|2010al ) , our geometry placed the 
inner edge of the wind at a relatively small inner radius and large outer radius, resulting in a 
comparatively large wind ope ning angle and ma king the mass distribution of the wind significantly 
different to that employed bv lSim et al.l (|2010al ). Our chose n geometry result ed in a relatively high 
ionization parameter compared to the model discussed by ISim et al.l ()2010al ) , owing to a reduced 
mean gas density and also less geometric dilution of the ionizing flux at the inner edge of the wind. 
Our chosen parameters also resulted in a faster flow. The higher ionization state of the gas results 
in a lower fraction of Fe xxv and Fe XXVI ions and so reduc es the dep t h of th e absorption features 
predicted by our solutions, compared to those predicted by ISim et al.l (|2010al ). The faster velocity 
of our flow also broadens the predicted absorption features, making them less promin ent in the 



model spectra. To achieve absorption features of a similar depth to those predicted by ISim et al 
(|2010al ). we would need to view the wind model presented here at higher angles than were found 



in the spectral fitting. In summary, the depth of Fe xxv and Fe XXVI absorption lines in the 
model spectra depend on the wind geometry, mass outflow rate and viewing angle, and there is no 
inconsistency between model predictions and the inclination angles found in our fits. 

Future observations with higher spectral resolution and signal-to- noise may be able to resolve 
the subtle spectral features in the model and allow all wind parameters to be determinable. In 
addition, soft band emission lines, such as O Ka, Si Ka and S Ka, may help to constrain some of 
the wind parameters. 



6. Conclusions 

We have taken a sample of sources that show moderately broad components of Fe K-shell 
emission, and fit them using a model for a Compton-thick outflow. All moderately broad Fe K 
emission components are well described by the wind model, demonstrating the general applicability 
of such models to the interpretation of X-ray emission lines. Such an approach complements studies 
of X-ray winds based on X-ray absorption features. Solutions found for the sample sources also 
show that 50 - 85 % of the flux in the 15-50 keV band has been transmitted through or scattered 
by the wind. 

We conclude that the X-ray signatures of a Compton-thick wind have a profound effect on the 
observed X-ray spectra of AGN, even when the nuclear system is not viewed through the wind. In 
principle, fits to X-ray spectra of AGN can provide important constraints on the wind mass loss 
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rate, geometry and orientation of the nuclear wind. Future work will extend this study to more 
complex sources, likely viewed at somewhat higher inclination angles. 
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Fig. 1. — The 2-10 keV fits: The black markers represent the data, the green curve represents 
the wind model component, and the red curve represent the total spectrum. Two observations are 
shown overlaid for Fairall 9: for clarity, the second observation (OBSID 7050603010) has the data, 
wind model, and total spectrum represented by cyan, orange, and magenta, respectively. 
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Fig. 2. — In the 0.6-50 keV fits: The black markers represent the data, the green curve represents 
the wind model component, and the red curve represent the total spectrum. Two observations are 
shown overlaid for Fairall 9: for clarity, the second observation (OBSID 7050603010) has the data, 
wind model, and total spectrum represented by cyan, orange, and magenta, respectively. 



Table 1. Parameters for Model Grid 



Parameter Values 



logio M w /M E dd 


-0.5 


-0.1 


logio R mi „/(100 Kg) 


-0.5 


0.0 


logio d/R m i n 


0.0 


0.5 


Rmax /Rmin 


3.0 


9.0 


Fe abundance/solar 


1.0 


3.0 


r 


1.8 


2.4 


fv 


1 




Rv 


Rmax 

l 




v 







M/Mq 


10 7 




k 


-1 






0.01 





Note. — The parameters with multiple values 
are the wind model grid points, while the param- 
eters with a single value are assumed in the wind 
model. See section 2.1 for details. 
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Table 2. Wind Fits 



Object 


r 


M„ 1 


2 


7V r 3 




kT 5 


T 6 


Rt 7 

1 l jrac 


Nh 8 


v 2 /do f 










2-10 keV 












Ark 120 


9 n9 +0.02 
z - uz -0.03 


u - o *-0.14 


63l^ 


99 o + 5.8 
zz - z -14.2 


85li' 5 


- 




1 00+ ' 48 
±.uu_ 72 


9.8l 8 ' 5 


1.05/124 


Fairall 9 (702) 9 


1 qt;+0.02 
i - yo -0.04 


co+0.85 


38l* 2 


25-21^3 


85l 4 ° 


- 




3.89l°f° 


3.2l 7 / 


1.23/250 


Fairall 9 (705) 9 


2 02+ ' 01 
z - u -0.03 


0.781-°, 


II 


25.8l? 4 2 8 


II 


— 


— 


45 +0.47 


II 


— 


MCG-02- 14-009 


1 84+0-08 


- 82 -0.49 


29l 2 ° 


3.2t°;! 




— 


— 


07+I.I9 
-2.99 


10.71 24 / 


1.02/124 


Mrk 335 9 


o i 1+0.02 


"S4+0- 16 


65l 2 


1 9 e;+ 3 - 5 
lz -°-9.5 




- 


- 


1 n^+o- 13 

l.UO_g g 2 


3.6l 6 / 


1.62/122 


NGC 7469 


1.80l° 04 


87+ 1 69 

U - S ' -0.54 


461JI 


o n+17.5 
d - y -2 .47 


00+47 


- 


- 


80+ 4 ' 20 


9.711.1 


0.76/124 


SWIFT J2127.4 


9 n9 +0.03 


0.461^° 


69l 2 


13 4+ 3 -4 
0.6-50 


keV 9 


- 


- 


40+ ' 36 




0.97/124 


Ark 120 


1 Q7+0.03 
l- y ' -0.03 


0.961°-? 


641^ 


17 7 +4.3 
L '-'-13.7 




< 11 


< 1.0 


0.7S±g;35 


10.2l«-8 


1.26/166 


Fairall 9 (702) 


1 qo+0.02 


u - oo -0.15 


65l 6 


23 1+ 4 ' 4 


85+° 


< 8.7 


< 2.5 




7 4+ 4 -6 


1.34/332 


Fairall 9 (705) 


9 n9 +0.02 
z - uz -0.03 


52+ ' 13 


II 


24.5l5 3 5 9 


11 


< 8.5 


< 2.3 


1.511^ 


11 




MCG-02- 14-009 


1 g + - 03 
1 -™-0 .02 


, 71 +0.«5 
1 - ' -0.82 


43±\° 


4 6+1-4 
^■"-1.9 


85l 3 " 






4.48l 3 ;f 7 


9.2l°- 5 


1.20/170 


Mrk 335 


^• uu -0.02 


0-5l±g;iS 


7ll 2 


1 q+0-8 


63l|| 


< 13.5 


< 2.0 


12+ - 08 


*-'-o.i 


1.88/166 


NGC 7469 


, sr -+0.02 

i - 8tl -o.oi 


0-86li;« 


44ll° 


9 O+0.8 
z -°-0.7 


85l 10 


< 13.9 


< 1.2 


46+ ' 24 
u -* D -0.28 




1.05/166 


SWIFT J2127.4 


9 n7 +0.02 
z - u '-0.02 


0.33l° 40 


69l 2 2 




85±° 






83+ ' 18 
u.oci_ 23 


79-Olij 


1.55/169 



Note. — Errors are calculated at 90% confidence. Values in italics represent ranges truncated by the limit imposed during the fit. 

1 The flow rate in M^;^, the Eddington mass-loss rate 

2 The observer's viewing angle through the wind, in degrees 

3 Normalization of the PEXMON component in units 10 — 3 photons cm _2 s _1 

4 The angle of inclination of the PEXMON component, in degrees 

5 COMPTT Plasma temperature, in keV 

6 COMPTT Plasma optical depth 

7 The calculated reflection scaling factor between PEXMON and the intrinsic powerlaw continuum 

8 Galactic column density in units of 10 20 atom cm~ 2 , with the values from the DL HI in the Galaxy survey as a lower limit 
9 Sec text for details 
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Table 3. 2-10 keV Luminosity 



Object Lofts 1 L x 2 L x /L Bdd 



Ark 120 3 


7.5E43 


2.1E44 


0.01 


Fairall 9 3 (702043010) 


1.2E44 


2.8E44 


0.03 


Fairall 9 3 (705063010) 


1.3E44 


2.6E44 


0.03 


MCG-02-14-009 4 


7.8E42 


1.6E43 




Mrk 335 5 


2.3E43 


5.9E43 


0.02 


NGC 7469 3 


1.3E43 


2.5E43 


0.03 


SWIFT J2127.4+5654 6 


1.6E43 


4.6E43 


0.03 



lr The observed 2-10 keV luminosity in erg s _1 
2 The intrinsic 2-10 keV luminosity in erg s _1 

3 Black hole mass estimate from I Woo &: Urrvl j2002h 

4 Black hole mass estimate not found in the literature 

5 Black hole mass estimate from lGrier et al". I l l2012l) 

6 Black hole mass estimate from lMalizia et"al I l l200g|) 



